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We improve the numerical kludge waveform model introduced in Huerta & Gair (2011) [il m 
two ways. We extend the equations of motion for spinning black hole binaries derived by Saijo 
et al. 01 using spin-orbit and spin-spin couplings taken from perturbative and post-Newtonian 
(PN) calculations at the highest order available. We also include first-order conservative self-force 
^SJ . corrections for spin-orbit and spin-spin couplings, which are derived by comparison to PN results. 

We generate the inspiral evolution using fluxes that include the most recent calculations of small 

■ body spin corrections, spin-spin and spin-orbit couplings and higher-order fits to solutions of the 
, Teukolsky equation. Using a simplified version of this model in HI , we found that small body spin 

t I . effects could be measured through gravitational wave observations from intermediate-mass ratio 

■ inspirals (IMRIs) with mass ratio rj J> 10^'^, when both binary components are rapidly rotating. In 
this paper we present results of Monte Carlo simulations of parameter estimation errors to study in 
detail how the spin of the small/big body affects parameter measurement using a variety of mass 

• and spin combinations for typical IMRI sources. We have found that for IMRI events involving a 

moderately rotating intermediate mass black hole (IMBH) of mass lO^M© and a rapidly rotating 
central supermassive black hole (SMBH) of mass W^Mq, gravitational wave observations made 
with LISA at a signal-to-noise ratio (SNR) of 1000 should be able to determine the inspiralling 
IMBH mass, the central SMBH mass, the SMBH spin magnitude, and the IMBH spin magnitude 
to within fractional errors of ~ 10~^, 10~^, 10~*, and 10~^, respectively. LISA should also be able 
Jh^- to determine the location of the source in the sky and the SMBH spin orientation to within ~ 10"^ 

steradians. Furthermore, we show that by including conservative corrections up to 2.5PN order, 
systematic errors no longer dominate over statistical errors. This shows that search templates that 

■ include small body spin effects in the equations of motion up to 2.5PN order should allow us to 
^ ■ perform accurate parameter extraction for IMRIs with typical SNR~ 1000. 
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I. INTRODUCTION 

Developing waveform templates that accurately model the inspiral of stellar mass compact objects (CO) into 
supermassive black holes (SMBHs) has enabled us to shed light on the scientific payoffs that would be obtained if we 
were able to detect and characterize a large number of these events. Current estimates suggest that if a space-based 
detector like LISA Q made a single detection of these extreme- mass-ratio inspirals (EMRIs), we could determine the 
intrinsic parameters of the system, i.e., the masses of the binary and the spin of the central SMBH, to a fractional 
accuracy of ~ 10~^, the orientation of the spin of the SMBH and the position of the source in the sky to a resolution 
of ~ 10"'^ steradians and the luminosity distance to an accuracy of ~ 10% [1, [![. 

Most existing EMRI models ignore the spin of the smaller object and so recent work on source modelling of EMRIs 
has tried to determine the impact of ignoring the small body spin effects on source detection and parameter estimation. 
Barack & Cutler estimated that over the lastyear of inspiral, the inclusion of small body spin effects would affect the 
accumulated orbital phase by a few radians [4?| . The reason being the dimensionless spin parameter of the inspiralling 
CO is multiplied by a factor of r/ in the equations of motion [2], and hence for EMRIs of mass ratio t] ~ 10^^, the 
contribution of the CO spin to the binary's dynamics is highly suppressed. Furthermore, Burko @ has shown that 
omitting local self-force corrections in the modeling of spinning COs that inspiral into supermassive Schwarzschild 
BHs is not needed to accurately model the orbital evolution, and that such omission will only decrease our ability to 
determine the spin rate of the companion. 

We have recently explored the effect of spin using a different approach. We developed a waveform model that 
includes small body spin corrections and determined the mass-ratio threshold 77 at which the small body spin effects 
become relevant for the dynamics. We found that the spin of the inspiralling CO could be accurately determined 
from gravitational wave (GW) observations only for 77 ^ 10""^. We also showed that, in this regime, small body spin 
corrections become particularly important when both components of the binary are rapidly rotating. 

A mass ratio of 77 ^ 10~^ can no longer be described as an EMRI, but represents an intermediate-mass-ratio inspiral 
(IMRI). This distinction is a modelling one. For small mass ratios, 77 <J 10^'* (EMRIs), the gravitational waveforms can 
be modelled accurately using black hole perturbation theory (BHPT), while for near equal- mass systems, rj ;> 10^^, 
the waveforms can be computed using post-Newtonian (PN) techniques. IMRIs lie between these two limits, in a 
regime in which the velocities of the binary components are too large for PN theory to apply, but the mass ratio is 
sufficiently large that perturbation theory cannot be used. 

Our findings in [l| suggested that search templates that aim to accurately measure the properties of IMRIs would 
also have to include small body spin effects. There are two different types of system that could be observed as 
IMRIs and both of these involve intermediate mass black holes (IMBHs) with mass M ~ 10^-10'' M©. An IMBH that 
inspirals into a SMBH in the centre of a galaxy could be detected by LISA if the SMBH has mass in the appropriate 
range ^ 10''-10^Mq. Similarly, if a stellar mass compact object, a neutron star or black hole, inspirals into an IMBH 
of mass ~ 10^ — the emitted GWs might be observed by ground-based detectors such as Advanced LIGO 

or the Einstein Telescope Such events could occur in the dense stellar environments of globular clusters that 

contain IMBHs. Both types of system constitute IMRIs, but in this paper we focus on the former category. 

For LISA IMRIs, an accurate determination of the spin of the inspiralling CO will provide information about the 
formation channels of IMBHs. Numerical studies suggest that high observed spins would imply that much of the mass 
of these objects was obtained by gas accretion from a disc, whereas moderate spins would favour a scenario in which 
the IMBH formed as a result of a major merger between comparable-mass black holes, and low spins would imply 
that the IMBH acquired its mass through a series of minor mergers with smaller objects pH. [l^. Therefore, it is 
important to determine if templates that include small body spin effects at the highest order available will allow us 
to measure the spin distribution of IMBHs through IMRI observations. In this paper we extend the results of [l[ by 
including higher order corrections in the waveform model and exploring a wider variety of possible systems. 

To perform this analysis, we will make use of an improved version of the "numerical kludge" waveform scheme. In 
this framework, one combines an accurate prescription for the generation of the inspiral trajectory with a pseudo-flat 
spacetime wave generation formula. Even though this construction is inconsistent, it has been shown that numerical 
kludge waveforms capture the main features of Teukolsky-based waveforms [T^. In our improved scheme, we extend 
the spinning black hole binary model derived by Saijo et al by: (i) generating the inspiral trajectory using fluxes of 
energy and angular momentum that include BHPT calculations for spin-spin and spin-orbit couplings at the highest 
order available, and (ii) by amending the orbital phase evolution with conservative self-force corrections for spin- 
orbit and spin-spin couplings using the same method that we used to compute perturbative conservative corrections 
in 

The analysis carried out in suggests that the spin distribution of massive Kerr BHs is heavily skewed toward 
rapid spins below redshift z ^ 5. This analysis suggests that gravitational radiation generated by IMRI sources 
involving supermassive Kerr BHs may include small body spin effects. Hence, it is worthwhile carrying a detailed 
study of the influence of higher-order spin corrections for parameter estimation of IMRIs. For completeness in our 
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analysis, we also estimate the systematic errors that arise when we omit higher-order corrections in the waveform 
template. We find that by including the conservative corrections up to 2.5PN order in search templates, it should be 
possible to reduce the magnitude of systematic errors so that they do not dominate over the parameter errors arising 
from instrumental noise, even for strong sources. 

This paper is organized as follows. In Section we describe our updated kludge waveform template that includes 
the most recent results from both perturbative and PN calculations for spin-spin and spin-orbit calculations in a 
physically consistent manner. In section IIIII we summarize the wave generation model used in our studies. In 
Section IIVI we describe the Fisher Matrix formalism used to estimate parameter estimation errors. We use this 
formalism in Section |V] to explore how the accuracy of parameter measurement depends on the spin of the small/big 
body using a representative sample of binary systems and a variety of spin configurations. In Section IVll we estimate 
the systematic errors that arise from omitting higher order corrections in the waveform template using the formalism 
introduced in 15]. We summarize our findings in Section [VIII 



II. IMPLEMENTATION OF HIGHER-ORDER SPIN EFFECTS IN THE EQUATIONS OF MOTION OF 

SPINNING BLACK HOLE BINARIES 

In the first paper of this series [l|, we showed how to extend the kludge waveform model introduced in [l^ by 
including small body spin corrections. In this Section, we briefiy review this scheme and show how we can improve 
this model by including perturbative and PN corrections for spin-spin and spin-orbit couplings at the highest order 
available. Thereafter, we will augment the search template with conservative self- force corrections which we derive 
by comparison to PN results. The aim of this analysis is to provide a search template that includes the most recent 
calculations, both from perturbative and PN analysis, to shed light on the capabilities of future low-frequency detectors 
to measure the spin of IMBHs that inspiral into SMBHs, and to find out whether including conservative self-force 
corrections for spin-spin and spin-orbit couplings at the highest order available will be enough to prevent model errors 
to dominate over statistical errors. 

The equations of motion we shall use to build the IMRI numerical kludge waveform model are those derived in [2] , 
which describe the orbital evolution of a spinning CO that inspirals into a spinning SMBH. We shall assume that the 
spin angular momentum Si of the inspiralling CO is aligned with the central supermassive Kerr BH spin S2 and the 
orbital angular momentum L^. In this configuration, the orbital evolution can be modeled using Eqs. (l)-(9) of Q. 

The first improvement to the IMRI waveform model introduced in [jy] affects the generation of the inspiral tra- 
jectory of the spinning CO. In this analysis, we will use the radiation fluxes derived by Gair & Glampedakis 16| . 
augmented with accurate BHPT results that include small body spin corrections at the highest order available 17| . 
The expressions for energy and angular momentum flux become 

73^ 25_ \ flY^^ ^ ( 44711 , 33 -2 , 71 ^ /^l^^ 




higher order Teukolsky fits > , 



61 19 \ flV''^ ( 44711 33 2 59 \ /l^^ 



higher order Teukolsky fits > , (1) 



where r = p/M, and p stands for the Boyer-Lindquist radial coordinate, s = s/M — rjx, where x stands for the 
dimensionless spin parameter of the inspiralling CO, j] ~ fJ'/M represents the mass-ratio, and q is the dimensionless 
spin parameter of the central SMBH. 

It is worth pointing out that to derive the conservative self-force corrections we will make use of the terms explicitly 
included in the expressions for E and Jz in Eq. ([T]). However, the waveforms will be generated including the higher- 
order Teukolsky fits derived in [l^ . 

Furthermore, as shown by Tanaka et al. [iTj . to evolve circular equatorial orbits for a spinning particle in a Kerr 
background we can use the "circular goes to circular" rule of [ll], i.e., [l3l 
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where d4>/dt = r2(r) stands for the orbital frequency. We should bear in mind that this assumption is consistent with 
the energy and angular momentum loss rates at linear order in the spin of the particle [l7j . 
These considerations then allow us to evolve the radial coordinate as follows 

d7' ■ dr • 

'^^d^^=dj:^- 

We will generate the inspiral trajectory, r{t), using the expression for the angular momentum quoted in ([T|) which 
includes higher-order spin corrections and higher-order fits toTeukolsky data. 

The search template we have developed so far includes radiative self-force corrections at the highest order avail- 
able. Hence, for completeness in the modeling, we shall now include conservative self-force corrections. Because the 
conservative piece of the self-force accumulates and affects the phasing of the waveform over time, we need to include 
this effect as it could lead to several cycles of phase discrepancy in the kludge model over the inspiral. We take into 
account this effect by rewriting the orbital frequency fl in the following form, 



+ (4) 

ycodcsic 



dt \ dt 



where the first term on the righthand side of Eq. Q is implicitly defined in Eq. and S^l is a frequency shift 
which depends on the instantaneous orbital parameters, and which is still unknown for spinning particles that inspiral 
into Kerr SMBHs. At present, the self-force program has been successful in computing gravitational self-force data 
for bound eccentric geodesies in Schwarzschild space-time, and in developing numerical algorithms to calculate the 
self-force on a scalar charge moving in bound eccentric equatorial geodesies in Kerr space-time [T9l - [23j . The extension 
of these numerical algorithms from the scalar to the gravitational case is the main goal of the self-force community in 
the foreseeable future. Hence, for the time being, the computation of gravitational self-force data for spinning particles 
in Kerr space-time is beyond the scope of the self-force program. However, we can still make progress by using the 
conservative self-force corrections that have been derived within the PN formalism. In order to be consistent with the 
expressions for the radiative part of the self- force, see Eq. ([T]), we will use conservative PN self- force corrections up 
to 2.5PN order. 

To implement PN conservative self-force calculations in the kludge waveform, we will modify the kludge orbital 
frequency and its first time derivative. We will then compare these asymptotic observables with their PN counterparts. 
This comparison will enable us to determine coordinates between both formalisms, and compute the conservative 
pieces. We shall start by expanding Eq. ([2]) as follows 



29^x(-) +o(r) ), (5) 



which takes the following form when we include conservative corrections 
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(6) 



To modify the first order time derivative of the orbital frequency we will use 
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Note that for consistency with Eq. we have used r at 2.5PN order. 

To compute the first time derivative of the orbital frequency within the khidge formahsm, f2, we assume that the 
time evolution of the radial coordinate is given by Eq. ([5]), and choose a gauge in which the 77^ terms in r{t) that are 
not proportional to x vanish. Under these assumptions, the first time derivative of ([6]) is given by 
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where, as before, s = rjx- To relate the kludge and PN formalisms we propose the following coordinate transformation 



bo + 62.1 ( ^ j + (63.2 + h.sq + xhA) 



3/2 



+ (64.1 + 64.2 gx) 



+ {b5.2 + b5.3q + b5,4X) 



I'M 



5/2 



(10) 



where R stands for the PN semi-major axis. 

For consistency with the analysis outlined above, we will now use PN expressions for the frequency and its first 
time derivative which include spin-spin and spin-orbit corrections at 2.5PN order [23 - [27l | 
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where m = M + /i, L is a, unit vector directed along the orbital momentum, X — Xi — S\/ fi^, Q — X2 — S-^/M"^, 
dm = fj, — M, and we have used the spin variables 
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so that Sf 



S ■ L,Y.f = T, ■ L. 



The first order time derivative of the orbital frequency, including spin-orbit effects at 1.5PN and 2.5PN beyond the 
dominant approximation, is given by (24-i27| l 
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(milpN)^^'^, and the secularly constant spin variables S'^, S'^ are given by S'^ = 



and 




where x ^ 

il^ , with Sl = S"" • L, SI = • L. The expressions for S^, are taken from Eqs. (7.2a) and (7.2b) 

of [25], under the restriction that the equations of motion used to build the kludge waveform model are valid when 
the spin of the inspiralling CO is perpendicular to the equatorial plane, and parallel to the momentum of the central 
object, i.e., for spinning circular-equatorial, non-precessing binaries. Note that the constant ^ in Eq. (|14l) has been 
included to guarantee that the PN framework and the perturbative approach coincide in the test mass particle limit 
?7 — >■ 0. The spin-spin parameter a is given by 



a = ^ (-247X1 • X2 + 721i • xiL ■ X2) ■ (15) 

In order to facilitate the comparison between kludge and PN results, we need to re- write expressions (ITT]) . (ITTl) in a 
way that takes the small mass-ratio limit by writing m — A/(l -I- 77). Finally, to obtain the conservative corrections in 
Eq. ([5]), we substitute relation ([TU|) into Eqs. ^ and © and solve simultaneously for all the coefficients by comparing 
the kludge expressions for ft and Ct with their PN counterparts, i.e., (fTTj) and respectively. We find that that 
the non-vanishing coefficients are given by 
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These are all the ingredients we require to explore the importance of including higher-order spin effects in the dynamics 
of compact binaries. To recap, the IMRI kludge waveform model we have developed is based on the equations of motion 
derived by Saijo et al [2|, which describe the dynamic evolution of spinning compact binaries. We generate the inspiral 
trajectory by evolving the geodesic parameters, see Eq. ([3]), using the most accurate expression currently available 
for the flux of angular momentum, which includes perturbative corrections for spin-spin and spin-orbit couplings, 
and higher-order Teukolsky fits for the inspiral of spinless particles. We have derived an improved expression for the 
azimuthal frequency which includes small body spin effects, and conservative self-force correction for higher-order spin 
effects. This asymptotic observable takes the following form when we implement all these corrections, see Eq. (jT]), 
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This expression includes higher-order spin effects up to 2.5PN order. This is the highest order currently available for 
perturbative calculations that include spin-spin and spin-orbit corrections. Hence, for consistency in this analysis, we 
used PN calculations at the same order to amend the azimuthal frequency. 

In the following sections we will explore the influence of higher-order spin effects on source detection and parameter 
estimation. We will accomplish this using two different approaches: 

(i) In Section|Vl we estimate the accuracy with which the spin of the inspiralling CO can be measured, and how the 
accuracy of parameter measurement depends on the spin of the small/big object. We do this using astrophysically 
relevant IMRI systems and a variety of spin configurations for the components of the binary. 

(ii) In Section I VII we assess the importance of the conservative self-force corrections for signal detection and 
parameter estimation by turning these parameters on and off in the model. 

III. WAVEFORM MODEL 

Following the numerical kludge waveform model developed in [Tsj . and its subsequent extension in [l|, we will 
construct the numerical kludge waveforms combining flat-spacetime wave equations with the true inspiral trajectory 
developed in Section |TT1 To apply the wave generation formula to the inspiral trajectory, we first identify the Boyer- 
Linquist coordinates of the source with flat-space polar coordinates. This scheme is inconsistent, but it is able to 
reproduce true gravitational waveforms in a large part of parameter space [l^ . 

The small perturbation hij (t) can be expressed in the transverse-traceless (TT) gauge as follows 

= [PtkPji - ^^PtjPk^j and = 'q^j - nSj, (18) 

where rjij stands for the flat Minkowski metric, D is the distance to the source, hi is a unit vector in the direction of 
propagation, and I^^ is the inertia tensor [IsJ . 

With regard to the LISA's response function used in these studies, we notice that since most of the SNR of the IMRIs 
considered in this work accumulates at frequencies f < lOmHz, it is adequate to use the low-frequency approximation 
to the detector response which was introduced in [29|, and which is summarized in Section IIIA of m. 

IV. NOISE INDUCED PARAMETER ERRORS 

To estimate the accuracies with which GW observations will be able to determine the system parameters, 0*, we 
use the Fisher Matrix formalism in the limit of high SNR. Given the expected sensitivity of LISA, this is appropriate 
for the systems we will consider in Section |Vl For large SNR, the covariance of the posterior probability distribution, 
(r~^)'^ , gives the expectation value of the errors A6'' 

{M'M^) = [T-^Y^ + 0(SNR)-\ (19) 
where the Fisher Information Matrix F*-' is given by 
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where is the true value of the parameters of the source. For white noise, i.e., S'„(/) = const., Parseval's theorem 
allows us to write the inner product of two given signals, (p | q), as 2S~^ J2a I^oo Paii)Qait)dt (J|. Hence, following 
Barack & Cutler we can define the "noise-weighted" waveform as follows 



This convention enable us to find the following approximate expression for the Fisher matrix 



Tab^^y^j daK{t)dM)At. (22) 

Using the same prescription presented in Section V(a) of Barack and Cutler [3], the total LISA noise function, Sh{f), 
has three components: instrumental noise, confusion noise from short-period galactic binaries, and confusion noise 
from extragalactic binaries. 



V. PARAMETER ESTIMATION ERROR RESULTS 



We shall now make use of the Fisher Matrix formalism introduced above to estimate the accuracy with which 
LISA will be able to determine the spin of IMBHs that inspiral into SMBHs. The IMRI kludge waveform model 
used in these studies includes all the higher-order spin effects described in Section |lll We present the parameter error 
estimates for two astrophysically relevant types of IMRI, assuming three different spin configurations: rapid spins, 
moderate spins, and low spins. 

It is worth pointing out that the parameter estimation results we present in the following Sections were obtained 
from Fisher Matrices which were nicely convergent over several orders of magnitude in the offset used to compute 
the numerical derivatives (see Eq. (PUI) ). We also checked that the inverses of these matrices had good convergent 
behavior. 

We will estimate the noise-induced errors by fixing the values of the intrinsic parameters of the source and running 
a Monte Carlo simulation over possible values of the extrinsic parameters. To estimate the SNR associated with each 
event we use the relation 



/■tLSO 

SNR2 = 2 / hl{t)dt, (23) 

and assume that the distance to the source is D — IGpc. We use these Monte Carlo results to obtain the SNR 
distribution of each type of source. This information will allow us to find a 'typical' SNR value to renormalise the 
parameter estimation results. We describe this method in detail in the following Section. 



A. Determination of typical SNRs 



To determine the appropriate SNR at which we quote the parameter estimation errors for each type of source, we 
have carried out a Monte Carlo simulation over values of the extrinsic parameters of the source, while fixing the value 
of the intrinsic parameters. The events we consider were uniformly distributed in comoving volume out to a redshift 
of z = 1, and were assumed to be observed over a one year time window at the detector. Under these assumptions, we 
computed the SNR for each event, then looked at the SNR distribution of the detected events and chose the "typical" 
reference SNR as the median of the SNR distribution. The events considered in this paper could be detected by LISA 
at redshifts z > 1. However, parameter determination for these events would involve the use of models that describe 
the formation channels of IMBHs and the inspiral rates -information which is currently uncertain. Hence, we shall 
consider events which take place only at redshifts z <^ 1. 

To ensure that the events are observed by LISA over one year, we choose the initial radius of inspiral, pq, such that 
at the end of the observation time window the inspiralling IMBH has reached the innermost stable circular orbit. We 
use this information as input data in the Monte Carlo simulations. We also fix the redshifted masses of the sources 
to the values of the two IMRI systems we will consider later, i.e., 5 x IO^Mq -|- IO^Mq and W'^Mq + lO^Mg. 
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Figure [T] shows the normahzed cumulative distribution functions for the SNRs computed in this way, for binary 
systems with central SMBHs of redshifted mass IO^A/q and spin parameters q — 0.9, 0.3, 0.1. The inspiralling IMBHs 
have specific spin parameters x = 0.9, 0.3, 0.1, and redshifted mass ^ — 5 x lO'^Af0. We find that the median SNR 
for sources with spin parameters (x = 0.3, 0.1) is approximately the same as that obtained for x = 0.9, so we only 
show the SNR distribution for this representative case for each type of source. This is true for all three spins of the 
central SMBH. 

In Figure [2] we present the normalized cumulative distribution function for binary systems with redshifted masses 
10^ Mq + lO^Af0, for the same spin configurations used in Figure [T] 




SNR SNR 

FIG. 1. Normalised cumulative distribution function for the signal-to-noise ratio of a cosmological population of binary systems 
with redshifted masses = 5 x W^Mq, M = IO^Mq. The inspiralling IMBHs were taken to have spin parameter x = 0.9, 
while the spin parameter of the central SMBH is q = 0.9 (top panel), q — 0.3 (bottom left) and q — 0.1 (bottom right). 



The results presented in Figure[l]show that the "typical" SNR of sources with redshifted masses 5 x IO^Mq + IO^Mq, 
and spin parameters {q, x) = (0.9, {0.9, 0.3, 0.1}); (g, x) = (0.3, {0.9, 0.3, 0.1}); and {q, x) = (0.1, {0.9, 0.3, 0.1}) are 
1000, 500, 500, respectively. Similarly, Figure [2] suggests that the median SNRs of the binaries with component 
masses IO^Mq + IO'^Mq and spin parameters (q, x) = (0.9, {0.9, 0.3, 0.1}); {q, x) = (0.3, {0.9, 0.3, 0.1}) and (g, x) = 
(0.1, {0.9, 0.3, 0.1}) are 1000, 750, 750, respectively. We will use these reference values to renormalise the parameter 
estimation errors presented in Tables HHIVIII so that they represent the typical parameter estimation errors for a 
cosmological population of sources at z ^ 1. For the top panel of Figure [TJ SNR= 1000 corresponds to the 60-th 
percentile of the distribution rather the median, but we have chosen this reference value as it is a convenient reference 
point which is easy to readjust to other SNRs. The results in Figures [T] and [5] could also be used to renormalise the 
results to a different SNR reference value if required. 
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FIG. 2. Normalized cumulative distribution function for the signal-to-noise ratio of a cosmological population of binary systems 
with BHs of redshifted mass ^ — W^Mq, M = IO^Mq. The panels show results for three different combinations of the spin 
magnitudes of the central and inspiralling BHs, (g, x): namely: top panel (0.9,0.9); bottom-left panel (0.3,0.9); bottom-right 
panel (0.1,0.9). 



B. Parameter estimation results 

In this Section we present parameter estimation results for the 11 parameters — see Table H] — that characterise 
these systems. The results obtained in [i| suggested that for IMRIs with mass ratio rj ^ 10~^, LISA observations 
might allow us to measure the spin of an inspiralling IMBH to a fractional accuracy of ~ 10% if observed with an 
SNR of 1000, provided both components of the binary were rapidly rotating, i.e, ioi q — x — 0.9. In this section, we 
extend that analysis and present results that suggest that GW observations could still provide information about the 
spins of IMBHs even if the spin of the inspiralling IMBH is not high. 

In Table im we can see that for IMRIs with mass components 5 x 1O'^M0-|-1O^M0, GW observations could determine 
the spin of IMBHs to a fractional accuracy of ~ 10% at a fixed SNR=1000 when the central SMBH is rapidly rotating 
{q = 0.9) and the spin of the IMBH is moderate to high {x — 0.9, 0.3). This is the first relevant result of this paper, 
and constitutes an important extension to the analysis originally presented in [l|, which did not include parameter 
estimation results for binaries with spin parameters {q = 0.9, x = 0.3). It is worth pointing out that this new finding 
does not crucially depend on the inclusion of higher-order spin effects in the waveform template. We reach this 
conclusion by comparing the results of Table V in [T] with those of Table [Tl] in this paper, which suggest that the 
inclusion of higher-order spin effects will have a minor impact on the accuracy with which GW observations may 
determine the source parameters of 5 x + IO^Mq IMRIs. We shall see below that the inclusion of higher-order 

spin corrections in search templates do become important for parameter estimation of lO^M© -I- IO^Mq IMRIs. 

Table IhI also shows that LISA observations could also constrain the spin of IMBHs with x ~ 0.1, but only to an 
accuracy ^ 25%. Tables ITTTl IIVI show the corresponding results for cases in which the binary has the same component 
masses as in Table HIl but the central SMBH has spin parameter q = 0.3, 0.1, respectively. These Tables confirm that 
we are unlikely to accurately measure the spin of an IMBH that inspirals into an SMBH with moderate or low spin. 
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In fi 


mass of inspiralling object 


InM 


mass of central SMBH 


q 


magnitude of (specific) spin angular momentum of SMBH 


X 


magnitude of (specific) spin angular momentum of inspiralling object 


Po 


Initial radius of inspiralling object's orbit 




Initial phase of inspiralling object's orbit 


es 


source sky colatitude in an ecliptic-based system 


<f)S 


source sky azimuth in an ecliptic-based system 


Ok 


direction of SMBH spin (colatitude) 




direction of SMBH spin (azimuth) 




distance to source 



TABLE I. This table describes the meaning of the parameters used in our model. The angles {6s, 4>s) and {9k, 4'k) are defined in a fixed 
ecliptic-based coordinate system. 



Model 


Distribution of logj^Q(AX) in error, AX, for parameter X — 


ln(m) 


ln{M) 


q 


X 


po 




es 


0s 


ex 




ln(D) 


q = 0.9 
X = 0.9 


Mean 


-3.14 


-2.97 


-4.37 


-1.08 


-1.67 


-1.55 


-2.18 


-2.13 


-1.85 


-1.75 


-1.92 


St. Dev. 


0.076 


0.94 


0.078 


0.071 


0.085 


0.438 


0.467 


0.345 


0.434 


0.530 


0.385 


L. Qt. 


-3.15 


-2.96 


-4.39 


-1.10 


-1.65 


-1.90 


-2.43 


-2.35 


-2.13 


-2.10 


-2.11 


Med. 


-3.13 


-2.95 


-4.37 


-1.06 


-1.65 


-1.71 


-2.02 


-2.01 


-1.89 


-1.80 


-1.95 


U. Qt. 


-3.06 


-2.82 


-4.28 


-1.01 


-1.53 


-1.46 


-1.81 


-1.79 


-1.62 


-1.58 


-1.73 


q = 0.9 
X = 0.3 


Mean 


-3.12 


-2.94 


-4.38 


-1.05 


-1.54 


-1.45 


-2.06 


-2.03 


-1.84 


-1.71 


-1.91 


St. Dev. 


0.059 


0.073 


0.053 


0.048 


0.074 


0.439 


0.427 


0.455 


0.381 


0.432 


0.297 


L. Qt. 


-3.14 


-2.95 


-4.39 


-1.06 


-1.63 


-1.89 


-2.34 


-2.24 


-2.14 


-2.03 


-2.10 


Med. 


-3.13 


-2.93 


-4.38 


-1.04 


-1.55 


-1.45 


-1.98 


-1.93 


-1.86 


-1.78 


-1.93 


U. Qt. 


-3.12 


-2.89 


-4.37 


-1.03 


-1.52 


-1.46 


-1.76 


-1.72 


-1.63 


-1.51 


-1.74 


q = 0.9 

x = o.i 


Mean 


-2.70 


-2.51 


-3.98 


-0.65 


-1.22 


-1.46 


-2.08 


-2.06 


-1.79 


-1.66 


-1.86 


St. Dev. 


0.124 


0.141 


0.151 


0.163 


0.221 


0.378 


0.449 


0.347 


0.391 


0.457 


0.495 


L. Qt. 


-2.85 


-2.65 


-4.16 


-0.82 


-1.37 


-1.88 


-2.43 


-2.35 


-2.11 


-1.96 


-2.06 


Med. 


-2.64 


-2.45 


-3.92 


-0.63 


-1.22 


-1.52 


-2.03 


-2.01 


-1.81 


-1.73 


-1.89 


U. Qt. 


-2.48 


-2.32 


-3.79 


-0.46 


-1.12 


-1.37 


-1.77 


-1.74 


-1.63 


-1.43 


-1.63 



TABLE II. Summary of results of the Monte Carlo simulation of Fisher Matrix errors for spinning BE systems with masses 
/u = 5 X 10^ Mq, M = 10® Mq and for a central SMBH with spin parameter q = 0.9. We show the mean, standard deviation, 
median and quartiles of the distribution of the logarithm to base ten of the error in each parameter. Results are given for the 
kludge model with conservative corrections to 2.5PN order and quoted at a fixed SNR= 1000. 
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Distribution of logj^o(^-^) in error, AX, for parameter X = 



Model 


ln(m) 


ln(A/) 


9 


X 


Po 


4>Q 


es 


0s 


Ok 


<t>K 


ln(D) 


q = 0.3 
X = 0.9 


Mean 


-3.10 


-2.90 


-2.70 


-0.29 


-1.59 


-1.52 


-1.95 


-1.87 


-1.80 


-1.67 


-1.80 


St. Dev. 


0.080 


0.084 


0.075 


0.064 


0.081 


0.375 


0.413 


0.463 


0.348 


0.402 


0.250 


L. Qt. 


-3.12 


-2.95 


-2.76 


-0.34 


-1.63 


-1.61 


-2.18 


-2.07 


-2.01 


-1.93 


-1.96 


Med. 


-3.09 


-2.92 


-2.73 


-0.30 


-1.59 


-1.51 


-1.86 


-1.78 


-1.79 


-1.66 


-1.83 


U. Qt. 


-3.00 


-2.82 


-2.64 


-0.23 


-1.51 


-1.43 


-1.64 


-1.58 


-1.58 


-1.42 


-1.66 


q = 0.3 
X = 0.3 


Mean 


-3.10 


-2.91 


-2.66 


-0.28 


-1.59 


-1.53 


-1.92 


-1.85 


-1.77 


-1.63 


-1.78 


St. Dev. 


0.073 


0.098 


0.057 


0.062 


0.087 


0.291 


0.397 


0.441 


0.344 


0.390 


0.247 


L. Qt. 


-3.13 


-2.94 


-2.73 


-0.32 


-1.58 


-1.62 


-2.14 


-2.07 


-1.98 


-1.88 


-1.93 


Med. 


-3.11 


-2.90 


-2.64 


-0.28 


-1.56 


-1.52 


-1.82 


-1.78 


-1.75 


-1.62 


-1.81 


U. Qt. 


-3.01 


-2.85 


-2.62 


-0.24 


-1.53 


-1.42 


-1.60 


-1.57 


-1.55 


-1.39 


-1.64 


q = 0.3 

x = o.i 


Mean 


-3.07 


-2.87 


-2.63 


-0.20 


-1.59 


-1.45 


-1.91 


-1.85 


-1.77 


-1.62 


-1.78 


St. Dev. 


0.081 


0.097 


0.093 


0.088 


0.102 


0.439 


0.421 


0.434 


0.349 


0.376 


0.245 


L. Qt. 


-3.11 


-2.92 


-2.72 


-0.25 


-1.61 


-1.89 


-2.21 


-2.10 


-1.99 


-1.84 


-1.93 


Med. 


-3.03 


-2.84 


-2.65 


-0.21 


-1.57 


-1.45 


-1.84 


-1.78 


-1.77 


-1.63 


-1.80 


U. Qt. 


-3.00 


-2.79 


-2.61 


-0.17 


-1.49 


-1.46 


-1.61 


-1.56 


-1.55 


-1.37 


-1.63 



TABLE III. As Table nil but for a central SMBH with spin parameter q = 0.3. Results are quoted at a fixed SNR of 500. 



Model 


Distribution of logj^Q(AX) in error, AX, for parameter X — 


ln(m) 


ln(A/) 


9 


X 


Po 


</>o 


e.s 




Ok 


<t>K 


ln(D) 


q = 0.1 
X = 0.9 


Mean 


-3.09 


-2.92 


-2.61 


-0.16 


-1.62 


-1.58 


-1.94 


-1.85 


-1.73 


-1.62 


-1.76 


St. Dev. 


0.075 


0.072 


0.065 


0.063 


0.073 


0.380 


0.416 


0.434 


0.396 


0.443 


0.306 


L. Qt. 


-3.16 


-2.99 


-2.67 


-0.20 


-1.68 


-1.80 


-2.24 


-2.15 


-1.96 


-1.86 


-1.96 


Med. 


-3.09 


-2.92 


-2.61 


-0.16 


-1.62 


-1.61 


-1.92 


-1.87 


-1.74 


-1.63 


-1.79 


U. Qt. 


-3.02 


-2.86 


-2.57 


-0.12 


-1.57 


-1.42 


-1.68 


-1.63 


-1.49 


-1.39 


-1.59 


q = 0.1 
X = 0.3 


Mean 


-3.09 


-2.92 


-2.61 


-0.16 


-1.62 


-1.58 


-1.94 


-1.84 


-1.71 


-1.64 


-1.76 


St. Dev. 


0.092 


0.078 


0.075 


0.068 


0.083 


0.376 


0.409 


0.425 


0.402 


0.385 


0.350 


L. Qt. 


-3.16 


-2.99 


-2.67 


-0.21 


-1.69 


-1.85 


-2.22 


-2.11 


-2.02 


-1.85 


-1.95 


Med. 


-3.08 


-2.92 


-2.61 


-0.16 


-1.61 


-1.60 


-1.93 


-1.86 


-1.71 


-1.63 


-1.79 


U. Qt. 


-3.03 


-2.85 


-2.56 


-0.12 


-1.54 


-1.49 


-1.69 


-1.62 


-1.50 


-1.45 


-1.52 


q = 0.1 

x = o.i 


Mean 


-3.09 


-2.92 


-2.61 


-0.16 


-1.62 


-1.59 


-1.95 


-1.84 


-1.69 


-1.64 


-1.75 


St. Dev. 


0.081 


0.079 


0.072 


0.079 


0.076 


0.311 


0.376 


0.493 


0.394 


0.392 


0.492 


L. Qt. 


-3.16 


-2.99 


-2.67 


-0.21 


-1.69 


-1.85 


-2.22 


-2.11 


-2.02 


-1.85 


-1.95 


Med. 


-3.10 


-2.92 


-2.62 


-0.17 


-1.61 


-1.64 


-1.95 


-1.85 


-1.72 


-1.63 


-1.79 


U. Qt. 


-3.02 


-2.86 


-2.56 


-0.12 


-1.54 


-1.49 


-1.69 


-1.62 


-1.50 


-1.45 


-1.52 



TABLE IV. As Table Ull but for a slowly rotating central BH with spin parameter q — 0.1. Results are quoted at a fixed SNR 
of 500. Note that for slowly rotating SMBHs, the spin magnitude of the inspiralling IMBH has no noticeable influence on the 
accuracy with which the source parameters can be determined. 
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Table IVl shows that for IMRIs with mass components W^Mq + lO^M© and whose central SMBH is rapidly rotating, 
GW observations may allow us to obtain an accurate measurement of the spin distribution of IMBHs. These results 
show that at a fixed SNR= 1000, the spin of IMBHs could be determined to a fractional accuracy better than 10% 
for the spin combinations {q = 0.9, {x — 0.9, 0.3, 0.1}). 

Table IVII shows that for IMRIs with the same component masses considered in Table |Vl but with moderately 
rotating central SMBHs, LISA observations may enable us to constrain the spin of rapidly rotating IMBHs to an 
accuracy better than ~ 25%, at a fixed SNR=750. If, on the other hand, the IMBH is either moderately or slowly 
rotating, we could at least be able to constrain the spin parameter to an accuracy ~ 30%, at a fixed SNR=750. We 
also note that these results for the spin of the IMBH are a factor of ^ 2 better than those shown in Table IIIIl 

Finally, Table IVlIl shows that for IMRIs with mass components 10*Mq + 10^ Mq, and a slowly rotating central 
SMBH, LISA observations will not allow an accurate extraction of the spin of IMBHs. The best we could do is to 
constrain this parameter to better than 40%, at a fixed SNR=: 750. 

To recap, the analysis of these results suggests that, (i) for IMRI sources with mass ratio rj ^ 10^^, observed at 
a typical SNR= 1000, GW observations may allow us to measure the spin of IMBHs with rapid and moderate spins 
when the central SMBH is rapidly rotating to an accuracy ~ 10%; (ii) for IMRI sources with mass ratios t] ~ 10~^, 
observed at a typical SNR= 1000, LISA observations will enable us to determine the spin of IMBHs whether they are 
rapidly, moderately or slowly rotating, as long as the central SMBH is rapidly rotating, to an accuracy better than 
~ 10%; (iii) for IMRI sources with rj ~ 10~^, but with moderately rotating SMBHs, LISA observations will enable 
us to constrain the spin of IMBHs to an accuracy better than 30%, at a fixed SNR= 750, even if the IMBH is slowly 
rotating. 

This improved analysis also confirms our results in [l[ , namely, the determination of the spin of the inspiralling 
object through GW observations is mostly determined by the effect that spin couplings have on the orbital evolution. 
Our initial analysis in [l[ indicated that the spin coupling can be best determined when the central SMBH has high 
spin. However, in light of these new results, we conclude that while we still require the central object to rotate rapidly, 
the companion does not need to have high spin. Indeed, we have found that when we observe IMRIs with mass ratios 
T] ~ 10~^ with sufficient SNR, small body spin effects become significant enough on the orbital evolution of IMRIs 
that we can get accurate information about the spin distribution of slowly rotating IMBHs. 

These results are astrophysically meaningful, since the studies presented in [lj| suggest that IMBHs that grow 
through a series of minor mergers are expected to have typical spins x 0.3. Hence, GW observations will not only 
shed light on the existence of IMBHs, but may also provide information about their spin distribution -an astrophysical 
property that is determined by their channel of formation. These findings would, in turn, be useful to refine models 
of BH growth. 



Distribution of logj^Q(AX ) in error, AX, for parameter X — 



Model 


ln(m) 


ln(A/) 


q 


X 


Po 


00 


e.s 




Ok 


<t>K 


ln(D) 


q = 0.9 
X = 0.9 


Mean 


-3.05 


-2.87 


-4.09 


-1.12 


-1.51 


-1.58 


-2.10 


-2.04 


-1.86 


-1.72 


-1.89 


St. Dev. 


0.039 


0.062 


0.045 


0.052 


0.063 


0.333 


0.454 


0.470 


0.382 


0.411 


0.287 


L. Qt. 


-3.06 


-2.88 


-4.11 


-1.15 


-1.52 


-1.80 


-2.36 


-2.30 


-2.11 


-1.98 


-2.05 


Med. 


-3.04 


-2.87 


-4.09 


-1.11 


-1.51 


-1.63 


-1.97 


-1.93 


-1.85 


-1.77 


-1.93 


U. Qt. 


-3.01 


-2.85 


-4.07 


-1.09 


-1.49 


-1.48 


-1.74 


-1.72 


-1.67 


-1.55 


-1.75 


q = 0.9 
X = 0.3 


Mean 


-3.03 


-2.87 


-4.06 


-1.04 


-1.48 


-1.58 


-2.07 


-1.99 


-1.87 


-1.73 


-1.89 


St. Dev. 


0.096 


0.104 


0.101 


0.097 


0.085 


0.311 


0.445 


0.494 


0.369 


0.399 


0.283 


L. Qt. 


-3.05 


-2.89 


-4.09 


-1.09 


-1.54 


-1.79 


-2.29 


-2.24 


-2.11 


-1.96 


-2.04 


Med. 


-3.04 


-2.88 


-4.10 


-1.09 


-1.49 


-1.60 


-1.98 


-1.91 


-1.87 


-1.71 


-1.92 


U. Qt. 


-3.02 


-2.85 


-4.01 


-0.98 


-1.38 


-1.45 


-1.72 


-1.69 


-1.63 


-1.48 


-1.73 


q = 0.9 

x = o.i 


Mean 


-2.92 


-2.79 


-4.00 


-0.96 


-1.39 


-1.50 


-2.07 


-2.00 


-1.88 


-1.73 


-1.89 


St. Dev. 


0.109 


0.134 


0.074 


0.059 


0.086 


0.314 


0.443 


0.501 


0.366 


0.403 


0.282 


L. Qt. 


-3.02 


-2.83 


-4.09 


-1.05 


-1.40 


-1.70 


-2.31 


-2.21 


-2.08 


-1.95 


-2.05 


Med. 


-2.91 


-2.76 


-4.01 


-1.01 


-1.38 


-1.53 


-1.96 


-1.89 


-1.88 


-1.71 


-1.92 


U. Qt. 


-2.84 


-2.71 


-3.98 


-0.99 


-1.34 


-1.41 


-1.71 


-1.67 


-1.65 


-1.49 


-1.73 



TABLE V. As Table HH but for an inspiralling BH with mass ^ = lO^M©. Results are quoted at a fixed SNR of 1000. 
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Distribution of logj^o(^-^) in error, AX, for parameter X = 



Model 


ln(m) 


In(Af) 


q 


X 


Po 


4>Q 


es 


</>s 


Ok 


<t>K 


ln(D) 


q = 0.3 
X = 0.9 


Mean 


-2.99 


-2.84 


-2.70 


-0.63 


-1.50 


-1.59 


-2.06 


-1.97 


-1.87 


-1.72 


-1.88 


St. Dev. 


0.050 


0.054 


0.039 


0.063 


0.055 


0.313 


0.430 


0.456 


0.358 


0.392 


0.274 


L. Qt. 


-3.07 


-2.90 


-2.72 


-0.67 


-1.51 


-1.78 


-2.32 


-2.22 


-2.08 


-1.96 


-2.03 


Med. 


-3.00 


-2.83 


-2.69 


-0.61 


-1.49 


-1.61 


-1.97 


-1.88 


-1.87 


-1.73 


-1.90 


U. Qt. 


-2.93 


-2.79 


-2.65 


-0.57 


-1.48 


-1.46 


-1.73 


-1.70 


-1.63 


-1.45 


-1.72 


q = 0.3 
X = 0.3 


Mean 


-2.99 


-2.82 


-2.70 


-0.55 


-1.49 


-1.60 


-2.05 


-1.97 


-1.86 


-1.72 


-1.88 


St. Dev. 


0.089 


0.095 


0.059 


0.062 


0.119 


0.317 


0.420 


0.440 


0.359 


0.392 


0.274 


L. Qt. 


-3.06 


-2.90 


-2.72 


-0.57 


-1.53 


-1.79 


-2.30 


-2.18 


-2.08 


-1.96 


-2.03 


Med. 


-3.00 


-2.83 


-2.70 


-0.54 


-1.48 


-1.56 


-1.95 


-1.87 


-1.86 


-1.72 


-1.92 


U. Qt. 


-2.92 


-2.76 


-2.68 


-0.52 


-1.47 


-1.44 


-1.73 


-1.68 


-1.61 


-1.45 


-1.72 


q = 0.3 

x = o.i 


Mean 


-2.88 


-2.69 


-2.68 


-0.50 


-1.49 


-1.59 


-2.06 


-1.97 


-1.87 


-1.73 


-1.87 


St. Dev. 


0.049 


0.053 


0.103 


0.099 


0.298 


0.354 


0.407 


0.456 


0.358 


0.393 


0.275 


L. Qt. 


-2.93 


-2.72 


-2.72 


-0.55 


-1.61 


-1.78 


-2.29 


-2.22 


-2.08 


-1.95 


-2.05 


Med. 


-2.87 


-2.70 


-2.70 


-0.51 


-1.53 


-1.63 


-1.97 


-1.89 


-1.87 


-1.72 


-1.90 


U. Qt. 


-2.84 


-2.67 


-2.67 


-0.49 


-1.43 


-1.46 


-1.71 


-1.66 


-1.63 


-1.45 


-1.72 



TABLE VI. As Table |Vl but for an inspiralling IMBH with spin parameter q — 0.3. Results are quoted at a fixed SNR of 750. 
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TABLE VII. As Table |Vl but for an inspiralling IMBH with spin parameter q = 0.1. Results are quoted at a fixed SNR of 
750. Note that as in Table ITVl when the central SMBH is slowly rotating, the determination of the source parameters has no 
strong dependance on the the spin magnitude of the inspiralling IMBH. 
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It is important to ask how much of an impact the additional terms we have included in the model have made 
to the parameter estimation results. The 'best case scenario', i.e., IMRIs with mass components IO^Mq + IO^Mq 
whose central SMBH is rapidly rotating, was not used for the original analysis reported in We can evaluate 
the importance of the higher-order spin corrections by repeating the calculation for this new case using the previous 
model described in 1] to obtain parameter error measurements for IMBHs with spins x — 0-9; 0-3, 0.1. The results 
are shown in Table IVIIII 
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TABLE VIII. As Table El but inchiding small body spin corrections up to 2PN order. Results are quoted at a fixed SNR of 
1000. 

A direct comparison between Tables fVl and IVllIl shows that the inclusion of higher-spin corrections in the waveform 
template increases our estimates of the accuracy with which GW observations might determine the parameters of 
these type of events. In particular, we notice that our ability to determine the spin of the inspiralling IMBH is 
improved by nearly a factor of two, as shown in Figure [3] 




0.2 0.4 0.6 0.8 

Spin 

FIG. 3. The plot shows the accuracy with which GW observations may determine the spin parameter of 10* IMBHs that 
inspiral into rapidly rotating IO^Mq SMBHs. We have used the median of a set of Monte Carlo simulations of Fisher Matrix 
errors to generate this plot. We assume that these sources are observed at a fixed SNR= 1000. The upper line has been 
obtained using a model with small body spin corrections up to 2PN order The bottom line includes all higher-order spin 
corrections, as described in this paper. 
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There is an additional reason to include higher-order spin corrections in the waveform template for this kind of 
sources. As we discussed in ilj, a waveform template that includes small body spin corrections up to 2PN order has 
systematic errors whose magnitude is roughly a factor ~ 4 larger than noise induced errors for intrinsic parameters. 
This initial analysis then suggested that in order to improve parameter error determination for IMRI sources, and 
enhance the reliability of these results, it was necessary to include higher-order spin corrections in the equations of 
motion to further reduce the ratio of systematic to noise induced errors. In the following Section we show that this 
approach is the right strategy to improve the accuracy of the waveform template. 



VI. MODEL-INDUCED PARAMETER ERRORS 



In the previous Section, we explored the capability of future low-frequency GW detectors, such as LISA, to determine 
the spin distribution of IMBHs through IMRI observations. We have shown that the inclusion of higher-order spin 
effects in the dynamics of compact binaries has an important effect, particularly for IMRIs with mass ratio ry ^ 10""^. 
In this Section, we explore the importance of including higher order PN corrections for parameter estimation. The 
strategy to perform this analysis consists of building templates that include all or only part of the corrections derived 
in Section HIl 

The IMRI kludge model in Eq. p7|) includes corrections up to 2.5PN order. In this Section, a search template that 
includes all these corrections will be taken as the 'true' waveform, /igr, whereas templates that include only some 
of the Ci coefRcients of Eq. ([6]), will be taken as approximate templates, /lAp. For instance, a template that includes 
corrections at 1.5PN order includes the coefRcients co,C2,C3 and C3.1 only. We explore the importance of including 
the corrections by changing the terms entering expression ()17p for the frequency, but leaving intact the radiative part 
of the waveform model described by Eqs. (H))-®- This approach is consistent with how the coefficients in Eq. p7)) 
were derived from comparison to PN models. An understanding of systematic errors is particularly relevant for IMRI 
sources, because, as shown in Section IVAl the typical SNRs range from a few hundred up to a thousand. Hence, it is 
possible that 'model' or systematic errors, which are SNR independent, could dominate the total parameter-estimation 
error. 

In order to assess whether systematic errors dominate, we will compute the ratio TZ of the systematic error to 
the noise-induced error. If this analysis indicates that 7?. ^ 1, then the statistical parameter error estimates that 
we presented in Section |V] are a reliable indicator of the precision to which we would measure source parameters in 
practice. If 7?, 3> 1 we would need to further improve the accuracy of the model, either by including higher-order 
PN corrections for non-spinning binaries or by including second-order conservative corrections in order to accurately 
measure parameters. To estimate the magnitude of the systematic errors, we will use the formalism introduced in 
p^ . At first order in the error A0*(n), we can identify two contributions. The first is due to noise in the detector, 
and the second can be related to the approximate nature of the waveform template, which introduces another source 
of error given by 



T-\d)y (a,hAp(0)|hGR(^)-hApW), 



(24) 



where {/igrI^*')} {/iapI^*')} represent the two manifolds in the vector spaces of data streams corresponding to 
the 'true' and approximate waveforms, respectively. Given the fact that: (i) we do not currently have 'true' waveform 
templates at hand and we can only experimentally determine the hAp(^) that is the best fit waveform for the data 
stream, s = hcpX^') + n; and (ii) we are unsure about the error = 9 — 9, since we do not know 9, we make progress 
by replacing, at leading order, hcpX^) — hAp(^) by hGR(6') — hAp{9) to obtain 



Ath^?'' - (r-\9)Y (^djhAp{9)\hGR{9) - hAp{9) 



(25) 



The Fisher matrix is to be evaluated using the approximate waveforms Tij{9) = {dih.Ap{9)\djhAp{9)). As discussed 
in [1^, this scheme works best when the waveform is rewritten using an amplitude-phase decomposition, i.e.. 



/i"(/)=A"(/)e^*°(/), 
where the amplitude A and phase ^' are given by Eqs. (56)-(58) in [l|. 



(26) 
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In this form, the first order approximation of Eq. (1251) becomes 



(27) 



at e 



In a previous study [l| , we found that IMRI waveforms that included conservative corrections up to 2PN order would 
have a ratio TZ of model to noise-induced errors of 72. ^ 4. This result suggested that it was necessary to improve the 
IMRI model to find out whether the inclusion of higher-order spin effects and their associated conservative corrections 
could further reduce this ratio. Here we show that by using conservative corrections up to 2.5PN order we are able to 
reduce the ratio TZ down to 7?. ~ 1 for the intrinsic parameters (mass, spin etc.). This would mean that the parameter 
estimation errors presented in the previous Section should be reliable. More importantly, this suggests that we will 
not need to go beyond 2.5PN to recover the source parameters accurately, even for strong IMRI sources. 
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TABLE IX. Summary of Monte Carlo simulation results for the ratio of model errors to noise-induced errors, for IMRIs 
including an IMBH of mass /i = 5 x 10^ M© and a central SMBH with spin parameter q — 0.9. We show the mean, standard 
deviation, median and quartiles of the distribution of the logarithm to base ten of the ratio for each parameter. Results are 
given for various comparisons, as indicated and described in the text. A comparison "A vs B" uses model A as the true 
waveform and model B as the search template. Note that the noise-induced errors are quoted at a fixed SNR= 1000. 
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Once the gravitational self-force results for inspirals in Kerr are available, we will be able to extend this model 
by using fully relativistic calculations. Nonetheless, we will still have to use perturbative and PN corrections for 
spin-spin and spin-orbit effects, since the scope of the current self-force program is focused on inspiraling objects that 
are non-spinning jl^-H^ . The results we present below shed light on the level of accuracy that a waveform template 
should have to enable accurate parameter reconstruction for loud GW sources with two spinning components. 

In Table ITXl and \X \ we present the ratio TZ of model errors to noise-induced errors for the same systems considered 
in Section fVl and for two different comparisons, 2.5PN vs 2PN, and 2PN vs 1.5PN. The results show that the ratio 
TZ becomes smaller as the approximate waveform, hAp, becomes closer to the "true" waveform, Hgr, and more 
importantly, that for the 2.5PN 2PN comparison 7?, ~ 1 for the intrinsic parameters. The error ratio is somewhat 
bigger for the extrinsic parameters, 7?, ~ 4. The fact that the error ratio has improved for the intrinsic parameters but 
not the extrinsic parameters is perhaps not surprising. Information about the intrinsic parameters comes primarily 
from the waveform phasing while that on the extrinsic parameters comes from amplitude modulations. The corrections 
we have added into the model in this paper are corrections to the phase evolution and therefore we would expect them 
to primarily impact errors in the intrinsic parameters, which is what we see. This suggests that we need to include 
conservative corrections up to 2.5PN order in IMRI search templates in order to recover accurate estimates of the 
intrinsic parameters, but that with templates of this order, model errors should no longer be a limiting factor. 

As discussed in Section|Vl IMRIs with mass components IO^Mq -|- 10^ Mq may provide the best chance to determine 
the spin distribution of IMBHs that fall into rapidly rotating SMBHs, even if the IMBH is slowly rotating. Table H 
shows that for these sources, the ratio TZ becomes significantly smaller for the 2.5PN — > 2PN comparison and that 
7?. ^ 1 for the intrinsic parameters as in the previous case, while again 7?. ^ 4 for the extrinsic parameter^. We 
expect to do better for systems of higher mass ratio, because the corrections that we are omitting in this comparison 



are proportional to the mass-ratio. Cutler & Vallisneri [15[ also presented results for systems with component masses 
IO^Mq -t- 10^ Mq, although their results were for binaries in which neither component had spin and they compared 
a consistent 3PN approximate waveform to a "true" 3.5PN waveform. They found that the error ratios were larger 
for intrinsic parameters TZ ^ 10 than for the extrinsic parameters TZ ^ a. few. The results we have here are in fact 
fairly consistent with their findings, in the sense that both the estimated error ratio for the extrinsic parameters and 
the magnitude of the systematic errors for the intrinsic parameters are very consistent between the two studies. The 
apparent discrepancy arises almost entirely from a difference in the magnitude of the statistical errors for the intrinsic 
parameters, which are one to two orders of magnitude smaller in their work. The inclusion of spin in the waveform 
model introduces correlations that tend to decrease the precision with which the system parameters can be measured 
and this is the explanation for the difference between the current work and that previous study. 

In summary, we have demonstrated that a waveform template that includes all currently available higher-order 
spin effects and conservative PN corrections for spin-spin and spin-orbit effects might be sufficient to enable us to 
accurately reconstruct the parameters of loud sources, since the relative importance of the 2.5PN — > 2PN change is 
negligible, even for the loudest IMRI sources considered here. 



Strictly speaking, 7?. ~ 1 for all the intrinsic parameters, except for the central SMBH spin parameter, for which 7?. < 2. 
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logj^o of the ratio TZ of model to noise-induced error for parameter X = 
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-0.289 


-0.418 


-0.318 


-0.369 


-0.315 


-0.299 


-0.305 


Med. 


0.050 


0.070 


0.458 


0.099 


0.061 


0.235 


0.355 


0.362 


0.516 


0.513 


0.502 


U. Qt. 


0.513 


0.492 


0.579 


0.550 


0.444 


0.918 


1.339 


1.095 


1.324 


1.337 


1.411 



TABLE X. As Table HXl but for IMRI sources with mass components IQ'^Mq + 10^ M©. Noise-induced errors are quoted at a 
fixed SNR= 1000. 
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VII. CONCLUSIONS 



This paper is the second in a series that sheds hght on the importance of including smah body spin correction in 
the modehng of extreme and intermediate mass ratio inspirals. In the first paper ^IJ, we showed that these effects are 
important for IMRI sources with mass ratios 77 > 10~^, and that future low- frequency GW detectors might provide 
an accurate measurement of the spin of IMBHs that inspiral into SMBHs, as long as both compact objects are rapidly 
rotating. In this paper we have improved our initial analysis in several important ways. First, we improved the 
waveform model by including the best information currently available from perturbative and PN results (23 - [27| 
to amend the equations of motion that describe the inspiral evolution of spinning compact objects into Kerr black 
holes We showed how to implement first order conservative corrections for spin-spin and spin-orbit effects at 
2.5PN order to compute the evolution of the inspiralling object's orbital frequency. Furthermore, we generated the 
trajectory of the inspiralling CO using the most accurate prescriptions available for the fluxes of energy and angular 
momentum which include higher-order fits to Teukolsky-based evolutions and perturbative results at 2.5PN order 
for spin-spin and spin-orbit couplings. 

We then used this improved model to refine our understanding of the capability of future low-frequency detectors 
to measure the spin distribution of IMBHs that inspiral into rapidly rotating {q — 0.9) SMBHs. Our first important 
finding was that for IMRI sources with mass components 5 x 10'^7\fQ + 10^ M©, GW observations at a fixed SNR=:1000 
could determine the spin of the IMBH to a fractional accuracy of ^ 10% for IMBHs with moderate to rapid spins. 
This is an important result in light of the current understanding on the formation channels of IMBHs. It is currently 
accepted that moderate spins would favour formation through mergers with comparable mass BHs, while rapid spins 
would favour growth through accretion Tllj . Hence, these results indicate that future low- frequency GW observations 
will provide valuable information on IMBH growth. 

Furthermore, we have shown that for IMRI sources with mass ratios ?; ^ 10^^, higher-order spin effects must be 
included in search templates for these events. This is because a search template that includes these corrections should 
allow us, through GW observations, to measure the spin of IMBHs that inspiral into rapidly rotating SMBHs to 
fractional accuracies better than ~ 10% at a fixed SNR= 1000. Our results also suggest that when the central SMBH 
is moderately rotating, GW observations would enable us to constrain the spin of inspiralling IMBHs to ^ 30%, 
at SNR= 750. Thus, events with 77 ^ 10~^ and rapidly rotating central SMBHs can be considered as "golden" 
IMRI sources, since they will provide a wealth of information on the existence, astrophysical properties, and physical 
processes that lead to the formation of IMBHs. 

The second part of our analysis was an important check on our results. We showed in Section IV Al that IMRI 
sources are intrinsically loud, which means that systematic errors, which are SNR independent, could dominate the 
total parameter error budget and limit our ability to obtain accurate parameter error estimates. If systematic errors 
were larger than noise-induced errors, our parameter errors estimates would not be trustworthy. In Section IVIl we 
estimated the magnitude of systematic errors that would arise when part of the higher order corrections in the search 
templates were omitted. These results showed that the relative error when using a 2PN template to detect a 2.5PN 
signal is small at fixed SNR, for all the systems considered in our analysis. Indeed, for the 2.5PN vs 2PN comparison, 
the ratio TZ of the model errors to the noise-induced errors is of order TZ <^ 1 for the intrinsic parameters, even for 
the loudest sources detected with SNR of 1000. This analysis indicated that by including higher-order spin effects, 
and their associated conservative corrections in search templates up to 2.5PN order would help us reduce the ratio 
of model to noise-induced errors for intrinsic parameters by a factor of ^ 4, as compared to search templates that 
include these corrections up to 2PN order [ll]. Hence, we conclude that to accurately reconstruct the parameters of 
strong IMRI sources, we might not need to go beyond 2.5PN order. This finding complements the study carried out 
in [l5| , where the authors showed that model errors could be a limiting factor in recovering the parameters of loud 
sources in which the components were non-spinning. In the current paper we have focussed on IMRI sources in which 
both components have spin and our waveform model is based on accurate equations of motion for circular-equatorial 
spinning BH binary systems, allowing it to be more readily applied in the strong-field. 

There is a second type of IMRI that could enhance our understanding of the astrophysical properties of IMBHs 
— LIGO IMRIs, in which a stellar-mass object (~ IM0) falls into an IMBH IOOOMq). These could be detected 
by future ground-based detectors, such as Advanced LIGO [13] or the Einstein Telescope [31]- Previous studies have 
explored the capability of the Einstein Telescope to determine the astrophysical properties of these sources when the 
inspiralling black hole is not spinning [l3|- In light of the findings presented in this paper, we see the need to 
improve this model by including small body spin effects, and recent results from fully accurate numerical relativity 
simulations of systems with mass ratios 77 ^ 10~^ [3^. The development of these search templates will improve the 
scientific payoffs that will be obtained through GW observations in the Advanced GW detector era. 
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